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Abstract 

We present a preliminary measurement of the time-dependent CP asymmetry for the neutral B- 
meson decay — > cpK^. We use a sample of approximately 227 million i?-meson pairs recorded 
at the T{4:S) resonance with the BABAR detector at the PEP-II i?-meson Factory at SLAC. We 
reconstruct the CP eigenstates (j)Kg and (pK^ where (p — > K^K^ , vr+vr^, and is observed 

via its hadronic interactions. The other B meson in the event is tagged as either a B^ or B^ from its 
decay products. The values of the CP- violation parameters derived from the combined (f)K^ dataset 
are S^k = +0.50 ± 0.25 (stat)lo;o4(syst) and C^k = 0.00 ± 0.23 (stat) ± 0.05 (syst). In addition, we 
measure the CP-violating charge asymmetry Acv{B~^ — > (j)K^) = 0.054 + 0.056 (stat) ±0.012 (syst). 
All results are preliminary. 
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1 INTRODUCTION 



Decays of B mesons into charmless hadronic final states with a (p meson are dominated by 6 ^ sss 
gluonic penguin amplitudes, possibly with smaller contributions from electroweak penguins, while 
other Standard Model (SM) amplitudes are strongly suppressed [T|. In the SM, CP violation arises 
from a single complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [2]. 
Neglecting CKM-suppressed contributions, the time-dependent CP-violating asymmetries in the 
decays (j)K^ and B^ — > J/tp are proportional to the same parameter sin 2/3 where 

the latter decay is dominated by tree diagrams. Since many scenarios of physics beyond the SM 
introduce additional diagrams with heavy particles in the penguin loops and new CP-violating 
phases, comparison of CP-violating observables with SM expectations is a sensitive probe for new 
physics. Measurements of sin2/3 in B decays to charmonium such as B^ J/ipK^ have been 
reported by the BABAR 4 and Belle collaborations, and the world average for sin2/3 is 0.731 it 
0.056 6^. In the decay B^ cpK^ the Belle collaboration measures sin2/3 = -0.96 ± CSOll];?? [ZJ, 
while the BABAR collaboration (with a sample of approximately 114 million BB pairs) measures 
sin2/5 = 0.47 ± 0.34(stat)+o:o6(syst) in the decays B° ^ (j)K^ and 5° ^ 0A'O. 

In the SM, neglecting CKM-suppressed contributions, the direct CP violation in B^ — > (l)K~^ [S], 
detected as an asymmetry Acv = {^^k- ~ ^ipK+)/(^(i>K- + ^if>K+) in the decay rates T^^i = 
T[B^ (j)K^), is expected to be zero; in the presence of large new-physics contributions to the 
b sss transition, it could be of order 1 ^O]- The BABAR collaboration measures (with a sample 
of approximately 89 milhon BB pairs) Acp{B^ (pK^) = 0.04 ± 0.09 ± 0.01 |TT] . 

In this paper we report preliminary measurements of the time-dependent CP asymmetry in 
the decay B^ — > (j)K^ and the charge asymmetry in the decay B+ (f)K^ based on a sample of 
approximately 227 million BB pairs collected at the T(45) resonance with the BABML detector |12j 
at the PEP-II asymmetric-energy e'^e~ storage ring located at the Stanford Linear Accelerator 
Center. 

2 THE BABAR DETECTOR 

The BABMl detector is described elsewhere ^|. The primary components used in the analysis 
are a charged-particle tracking system consisting of a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH) surrounded by a 1.5-T solenoidal magnet with an instrumented 
flux return (IFR), an electromagnetic calorimeter (EMC) comprised of 6580 CsI(Tl) crystals, and 
a detector of internally reflected Cherenkov light (DIRC) providing excellent charged K and vr 
identification ;14] in the momentum range relevant for this analysis. 

3 ANALYSIS METHOD 

From a B^B^ meson pair we fully reconstruct one meson, Bcp, in the final state (pK^ , and partially 
reconstruct the recoil B meson, Btag. We examine Btag for evidence that it decayed either as B^ 
or B^ (flavor tag). The asymmetric beam configuration in the laboratory frame provides a nominal 
boost of (3"f = 0.56 to the T(4S'), which allows the determination of the proper decay-time difference 
At = tcp — ttag using the vertex separation of the two neutral B mesons along the beam [z) axis. 
The decay rate f+(f-) when the tagging meson is a B^{B^) is given by 

e-l^*l/^BO 

f±( At) = — [iTVfS^K sin (AmrfAt) ^ C^k cos (Am^At)] , (1) 
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where r^o is the neutral B meson mean hfetime, Arrid is the B^-B^ osciUation frequency, and the 
CP eigenvalue is rjf = —1 (+1) for (pK^ {(pK^). The time-dependent CP-violating asymmetry is 
defined as Acp = (f+ — f-)/(f++f-). In the SM, decays that proceed purely via the b sss penguin 
transitions have CP parameters St^x = sin 2/3 and C(j,K = 0, where (3 = arg [—Vcd^cb/^td^tbl- Here 
Vik is the CKM matrix element for quarks i and k. 

4 EVENT RECONSTRUCTION 

The Bcp candidate is reconstructed in the decay mode cpK^ with cj) — > K~^K~; the is either a 
or a Kg decaying into tt^tt". We combine pairs of oppositely charged tracks extrapolated to a 
common vertex to form (p and Kg candidates. For the charged tracks from the (p decay we require 
at least 12 measured drift-chamber coordinates and a minimal transverse momentum of 0.1 GeV/c. 
The tracks must also originate from within 1.5 cm of the nominal beam spot in the plane transverse 
to the beam axis and ±10 cm along the z-axis. Tracks with momentum less than 0.7 GeV/c that 
are used to reconstruct the (j) nieson are distinguished from pions and protons via a requirement 
on the likelihood that combines dE/dx information from the SVT and the DCH. For tracks with 
higher momentum, dE/dx in the DCH and the Cherenkov angle and the number of photons as 
measured by the DIRC are used in the likelihood. The two-kaon invariant mass must be within 
15 MeV/c^ of the known mass 0. 

For tracks corresponding to Kg and -Btag daughters our requirements are less restrictive. A 
Kg vr^vr" candidate is accepted if its two-pion invariant mass is within 15 MeV/c^ of the known 
K^ mass W, its reconstructed decay vertex is separated from the (p decay vertex by at least 3 
standard deviations, and the projected angle between the line connecting the (p and K^ decay 
vertices and the Kg momentum direction, in the plane perpendicular to the beam axis, is less than 
45 mrad. 

We identify a -fC^ candidate like in our B^ — > J/ tpK'l analysis T3' either as a cluster of energy 
deposited in the electromagnetic calorimeter or as a cluster of hits in two or more layers of the 
instrumented flux return that cannot be associated with any charged track in the event. The K^ 
energy is not well measured. Therefore, we determine the K^ laboratory momentum from its flight 
direction as measured from the EMC or IFR cluster, and the constraint that the invariant cpK^ mass 
agrees with the known B^ mass. In those cases where the K^ is detected in both the IFR and EMC 
we use the angular information from the EMC, because it has higher precision. In order to reduce 
background from vr*^ decays, we reject an EMC K^ candidate cluster if it forms an invariant mass 
between 100 and 150 MeV/c^ with any other cluster in the event under the 77 hypothesis, or if it 
has energy greater than 1 GeV and contains two shower maxima consistent with two photons from 
a tt" decay. The remaining background of K^ candidates due to photons and overlapping showers is 
further reduced with the use of a neural network constructed from cluster shape variables, trained 
on Monte Carlo (MC) simulated B^ (p^L ^^'^ measured radiative Bhabha events, and tested on 
measured e'^e~ — > (p{^ ^s^l)7 B^ J/tpK^ events. 

5 EVENT VARIABLES 

The results are extracted from an extended unbinned maximum likelihood fit for which we param- 
eterize the distributions of several kinematic and topological variables for signal and background 
events in terms of probability density functions (PDFs) QS]. The selection keeps loose requirements 
in those variables to include ranges dominated by background, too. The background B candidates 
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come primarily from random combinations of tracks produced in events of the type e~^e~ — > qq, 
where q = u,d,s,c (continuum). Background from other B decay final states with and without 
charm is estimated with MC simulations. Opposite-CP contributions from the K~^K~K^ final 
state {K^K~ S-wave) are estimated with data using a moment analysis method to be less 
than 6.6% at a 95% confidence level and are treated as a systematic error. The shapes of event 
variable distributions are obtained from signal and background MC samples and high statistics 
data control samples. In many cases parameters describing these distributions are varied in the 
likelihood fit. 

Each Bcp candidate is characterized by the energy difference AE = E*^ — i^y/s and, except 

for (pK^, the beam-energy-substituted mass uies = ^J {^s + po ■ pb)"^ /Eq — jEj. The 

subscripts and B refer to the initial T(45) and the Bqp candidate, respectively, and the asterisk 
denotes the T{AS) rest frame. For signal events, AE' is expected to peak at zero and mss at 
the known B mass. We require AE < 0.08 GeV for 5° (pK^, and |AE| < 0.1 GeV and 
iTT-ES > 5.21 GeV/ for B^ — > (j)Kg. The (/)-meson signal in the KK invariant mass, rnxK, is 
described with a relativistic P-wave Breit-Wigner function with parameters obtained from data. In 
the fit we also use the helicity angle Oh, which is defined as the angle between the directions of the 
and the parent Bqp in the K~^K~ rest frame. The cos 9h distribution for pseudoscalar-vector 
B decay modes is cos^ Oh, and for the combinatorial background it is nearly uniform. 

In continuum events, particles appear mostly in two jets. This topology can be characterized 
with several variables computed in the T(4S') frame. One such quantity is the angle 9t between the 
thrust axis of the Bqp candidate and the thrust axis formed from the other charged and neutral 
particles in the event. We also use the angle Op between the Bcp momentum and the beam axis, 
and the sum of the momenta pi of the other charged and neutral particles in the event weighted 
by the Legendre polynomials Lq^Oi) and L2{9i) where Oi is the angle between the momentum of 
particle i and the thrust axis of the Bcp candidate. For B^ — > 4>K^ candidates, we combine these 
variables into a Fisher discriminant J- JHI- In this mode, background from other B decays is 
negligible, as demonstrated in MC simulation studies. 

More stringent criteria must be applied to suppress backgrounds in the case of B^ — > 
candidates, and we require | cos^tI < 0.8 and | cosOb\ < 0.85. We define the missing momentum 
Pmiss, calculated in the laboratory frame from the sum of beam momenta and all tracks and EMC 
clusters, excluding the candidate. We require the polar angle Omiss of the missing momentum 
with respect to the beam direction to be greater than 0.3 rad. The cosine of the angle between 
Pmiss and the direction. Ok, must satisfy cos Ok > 0.6. In the plane transverse to the beam 
direction, the difference between the missing momentum projected along the direction and the 
calculated momentum must be greater than —0.75 GeV/c. In the Fisher discriminant we replace 
I cos 9b\ by the cosine of the angle between the missing momentum and the from the (p decay. 
The dominant CP contamination is the mode B (pK*^, where the K*^ decays to K^i:^. In the 
likelihood fit we explicitly parameterize backgrounds from both charm and charmless B decays, 
differently for neutral and charged B mesons, as derived from MC simulations. 

All the other tracks and clusters that are not associated with the reconstructed B^ — > (j)K^ decay 
are used to form the -Btag, and its fiavor is determined with a multivariate tagging algorithm |19j . 
The tagging efficiency and mistag probability Wi in six hierarchical and mutually exclusive 
categories are measured from fully reconstructed B^ decays into the D^*^~X~^ {X~^ = vr"*", p"*", a^) 
and J/tp K*^ {K*^ — > K~^tt~) fiavor eigenstates (-Bflav sample). The analyzing power X]^=i^«(l ~ 
2wi)2 is (30.5 ± 0.4)%. 

A detailed description of the At reconstruction algorithm is given in Ref. ^Hl- The Bqp vertex 
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resolution is determined by the (p vertex. The average Az resolution is 190 /xm and is dominated 
by the tagging vertex in the event. Thus, we can characterize the resolution with the much larger 
^flav sample, which we fit simultaneously with the CP samples. The amplitudes for the Bcp 
asymmetries and for the -Bflav flavor oscillations are reduced by the same factor due to wrong-flavor 
tags. Both distributions are convoluted with a common At resolution function. Backgrounds are 
accounted for by adding terms to the likelihood, incorporated with different assumptions about 
their At evolution and resolution function |15j . 



6 MAXIMUM LIKELIHOOD FIT 

Since we measure the correlations among the observables to be small in the data samples entering 
the fit (the largest one is 13% between niES and for the signal, all others are below 7%), we 
take the probability density function Vj^ for each event j to be a product of the PDFs for the 
separate observables. For each event hypothesis i (signal, backgrounds) and tagging category c, for 
the (/.i^-o mode we define P/'^ = ViimEs) ■ Vi{AE) ■ Vi{T) ■ Vi{mKK) ■ ViicosOH) ■ Pi(At; cja*, c), 
for the (pKl mode = ViiAE) ■ Vi{T) ■ Vi{mKK) ■ Vi{cos 0h) ■ ViiAt; a/^t, c), and for the flavor 

sample Vj^ = ViiruEs) - Vi^At; a/\t, c). The a^t is the error on At for a given event. The likelihood 
function for each decay is then 



£=nexp -j^iv.e n 



(2) 



where Ni^c is the yield of events of hypothesis i determined by the fit in category c, and Nc is the 
number of category c events in the sample. The total sample consists of 135,315 Sflav) 4300 (pK^ 
and 8238 4>K^ candidates. The reconstruction efficiency for the (j)Kg mode is about 40% and and 
for the (pKl mode about 20%. From the fit we find 114 ± 12 (pK^ and 98 ± 18 4>K^ signal events. 
The signal yields in both the (pK^ channels agree well with our determination of the branching 
fraction for (pK^ Figure H shows the mss (Ai?) distribution for (pK^ (0-^l) events 

together with the result from the fit after applying a requirement on the ratio of signal likelihood 
to the signal-plus-background likelihood (computed without the variable plotted) to reduce the 
background. 

We determine the CP parameters and C^k along with 83 other unconstrained parameters: 
event yields in signal and background (18 parameters), distributions of kinematic and topological 
variables for signal and background (12), the signal efficiency per tagging category (6), the average 
mistag fraction and the difference between and B^ mistags for each tagging category in the 
signal (12), and the signal At resolution (17). The At parameters for the charmless B background 
are the same as for the cpK^ signal. For the B decays into charm final states we parameterize 
the At resolution (3) and the mistag fractions (12). Their parameters are shared with the -Bflav 
sample. The At resolution for the continuum background (3) is kept unconstrained in the (pK^ 
datasets. We fix r^o and Am^ to the world averages jBj. The determination of the mistag fractions 
and At-resolution parameters is dominated by the large -Bflav sample. The fit was tested with 
both a parameterized simulation of a large number of data-sized experiments and a full detector 
simulation. The likelihood of our data fit agrees with the likelihoods from fits to the simulated 
data. The fit was also verified with our J /ipK^ and J/ipK^ data samples. 

As a cross check the analysis was also performed using different selection criteria which we 
describe in turn. The invariant K^K'^ mass is required to be within 10 MeV/c^ of the known mass 
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Figure 1: Distribution of the event variable (a) tues for the (l)Kg final state and (b) AE for the 
cl)K^ final state after reconstruction and a requirement on the ratio of signal likelihood to the signal- 
plus-background likelihood, calculated without the plotted variable. The signal efficiency for the 
selection and likelihood requirements is 32% for (a) and 9% for (b). The solid line represents the 
fit result for the total event yield and the dotted line for the total background. The dash-dotted 
(lower) line in (b) represents the continuum background only. 



of the 4> meson and is not used in the likelihood fit. The Kg flight requirements are tightened. The 
same four-category multivariate tagging algorithm as was used for the earlier published analysis 
is used. Instead of the Fisher discriminant a multivariate algorithm j20j for continuum background 
suppression is used, which in the (pKg final state combines the same four variables. In the case of 
the ingredients are Lq, L2, Pmiss, cos 6b, and cos 6t- The algorithm is trained in the same 
way as the Fisher discriminant and tested on data control samples. The central values of S^x and 
C(j)K for both the cross-check analysis and the primary analysis were hidden until the analyses were 
complete. We measure values for S^k and C^pK in very close agreement for the cpKg and the cp^L 
sample, separately, and for the combined samples. 

In the measurement of the CP-violating charge asymmetry in the decay 4>K~^ the se- 

lection of the (f) meson candidate is identical. For the candidate from the decay the track 
requirements are the same as for the (p daughters but we apply a more restrictive kaon identification 
criterion. We use the same set of event variables as for the ipK^ channel. The likelihood is the 
same as in Eq. I^with c corresponding to the two charge categories in signal and continuum back- 
ground. The total sample consists of 6654 (pK~^ candidates and from the fit we find 400 it 23 signal 
candidates. Figure [21 shows the distribution of niEs and rrixK with the result of the likelihood fit 
superimposed. We do not observe a significant asymmetry in Monte Carlo or in the continuum 
background data. 
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Figure 2: Distribution of the event variable (a) mES and (b) rnxK in the (j)K'^ final state after 
reconstruction and a requirement on the likelihood calculated without the plotted variable. The 
efficiency for the selection and likelihood requirements is 37% for (a) and 40% for (b). The solid 
line represents the fit result for the total event yield and the dotted line for the background. 



7 SYSTEMATIC STUDIES 

We consider systematic uncertainties in the CP coefficients and C^k due to contributions from 
final states with opposite CP (+0.06 for S^pxi ±0.02 for C^k)-, the parameterization of PDFs 
for the event yield in signal and background (±0.01, ±0.01), CP asymmetry of the background 
(±0.02, ±0.01), the assumed parameterization of the resolution function (±0.02, ±0.01), a 
possible difference in the efficiency for and B^ (±0.01, ±0.02), the fixed values for Am^ and tb 
(±0.00, ±0.01), the beam-spot position (±0.01, ±0.01), and uncertainties in the SVT alignment 
(±0.01, ±0.01). The bias in the coefficients due to the fit procedure is included as uncertainty 
(±0.01, ±0.01) without making corrections to the final results. We estimate errors due to the 
effect of doubly CKM-suppressed decays jU to be (±0.01, ±0.03). We add these contributions in 
quadrature to obtain the total systematic uncertainty. 

For the measurement of the charge asymmetry Acv we estimate the uncertainty due to charge 
asymmetries in tracking and particle identification to be 0.011. We also consider the systematic 
error due to uncertainties in the parameterization of the signal Fisher PDF (0.005) and B back- 
ground content (0.002). We add these contributions in quadrature to obtain the total systematic 
uncertainty. 

8 RESULTS 

The simultaneous fit to the (f)K^ and flavor decay modes yields the preliminary result 

S^K = ±0.50 ± 0.25 (stat)+°;°J(syst), 
C^K = 0.00 ±0.23 (stat) ±0.05 (syst). 
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The preliminary results in the channel — > (j)Kg alone are S^k = 0.29 ± 0.31 and C^j^ = 
-0.07 ± 0.27, and in the channel B° (l)K% S^k = 1-05 ± 0.51 and C^k = 0.31 ± 0.49, with 
statistical errors only. Figure |21 shows the At distributions of the B^- and the S'^-tagged subsets 
together with the raw asymmetry, for (fiK^ and (f^K^ events separately, with the result of the 
combined time-dependent CP-asymmetry fit superimposed. 

The preliminary value of the charge asymmetry in B^ — > (l)K~^ is 

Acv = 0.054 ± 0.056 (stat) ± 0.012 (syst). 

9 CONCLUSION 

In the decay B^ (j)K^ we measure preliminary values for S(f,K and C^k in the time-dependent 
CP asymmetry that are in close agreement with our previously published values 8 . Our value of 
S^K agrees within one standard deviation with the value of sin2/? in the B^ {cc)K^ decays |19) . 
We do not observe a significant charge asymmetry in the mode B^ cpK^ . 
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each final state, a requirement is applied on the event likelihood to suppress background. 
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